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We perform a combined analysis of antikaon-nucleon scattering cross sections and the recent SIDDHARTA kaonic hydrogen data 
in the framework of a coupled-channel Bethe-Salpeter approach at next-to-leading order in the chiral expansion of the effective 
potential. We find a precise description of the antikaon-proton scattering amplitudes and are able to extract accurate values of the 
scattering lengths, 
■A(1405). 
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1. Introduction and summary 

■ With the recent precise measurement of the characteristics of 
[kaonic hydrogen by the SIDDHARTA collaboration [H, an ac- 
curate determination of the so important antikaon-nucleon scat- 
^tering amplitude is now possible. The appropriate framework to 
■perform this task is unitarized chiral perturbation theory, which 
"combines the strictures from the chiral SU(3) dynamics of QCD 
■with coupled channel effects, that e.g. generate the much dis- 
[cussed A(1405) resonance, as first pointed out by Dalitz and 
■Tuan From earlier studies by various groups, it is already 
"known that simply taking the leading order chiral interactions 
■in the effective potential of the respective scattering equation 
is insufficient to achieve the desired accurate theoretical de- 
scription, see e.g. Refs. 13,11,11]. In fact, Ikeda et al. 
'have performed such a combined analysis based on the next-to- 
leading order chiral effective meson-baryon Lagrangian, nicely 
demonstrating that indeed a more precise description of the 
'K p and K n interaction arises. Here, we perform a similar 
[analysis, but in contrast to Refs. (Sill], we use a Bethe-Salpeter 
•framework without an on-shell approximation for the interme- 
diate meson-baryon states (as described in more detail below). 
The framework we use has already been successfully applied to 
pion-nucleon scattering in the s-waves ]8] and thus it is evident 
to extend this analysis to antikaon-nucleon scattering. We also 
point out that the constraints from SIDDHARTA on the kaon- 
deuteron scattering length have been investigated in ]|9|]. 

The main results of our investigation can be summarized as 
follows: 

• Fitting the scattering data for K p — » K p, SS'n, S^tt"^, 
and IPt^ for laboratory momenta piab < 300 MeV together 
with the SIDDHARTA data allows for a good descrip- 
tion of the antikaon-proton cross section data (cf. Fig. [U 
and an accurate determination of the scattering lengths, cf. 
Eq. ©. 



We can give a precise prediction for the real and imaginary 
part of the K p — » K p scattering amplitude for center-of- 
mass energies 1330 MeV < Wcms < 1450 MeV, cf. Fig.S 

We have investigated the two-pole structure of the A( 1405) 
1 in . While the first pole is in agreement with other 
determinations, we find the real part of the second pole at 
larger energies than usually obtained. We trace this back to 
the fitting procedures in other works that restrict the NLO 
contributions to come out close to the solution given by the 
Weinberg-Tomozawa term, while we do not impose such 
a restriction. 



2. Framework 

The starting point of our work is the recent analysis of nN 
scattering, presented in Isl]. In this section we describe the basic 
ingredients of this approach. 

In chiral perturbation theory the meson-baryon interaction at 
the leading chiral order is encoded in the following Lagrangian 



D/F 



where (. . .) denotes the trace in flavor space, D^Z? 
j[[u^, dfju], B], mo is the baryon octet mass in the chiral SU(3) 
limit, and D, F are the axial coupling constants. The relevant 
degrees of freedom are the Goldstone bosons described by the 
traceless meson matrix 0, which is included in the above La- 
grangian via :- exp{i(p/F{)) and m'^ := iu^ d^u- iu&^u^ . Here, 
F() is the meson decay constant in the chiral limit. The low- 
lying baryons are collected in the traceless matrix B. We set the 
external currents to zero except for the scalar one, which is set 
equal to the quark mass matrix, i.e. s - M:- diag(OT„, mj, m j. 
We furthermore use;\f± := u^x'^^ ± m^'m ^Vidx '■- ^Bqs, where 
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the constant Bq is related to the quark condensate in the chiral 
limit. 

Starting from the covariant derivative D^B, the so-called 
Weinberg-Tomozawa term can be derived. This term domi- 
nates the s-wave interaction near the thresholds, therefore in 
most chiral unitary approaches the meson-baryon interaction 
is restricted to this term. Secondly, a meson can couple to a 
baryon via the axial vector current ~ D, F, generating the s- 
and M-channel exchanges of the intermediate baryons. The in- 
clusion of these so-called Born graphs in the driving term of 
the Bethe-Salpeter equation leads to conceptional and practical 
difficulties, which are described in detail in Ref. i8l ll2tl . The 
latter are usually overcome, making use of the on-shell approx- 
imation or via projection of the kernel to the s-wave, see e.g. 
Ref. and Ref. |6] for a more recent study. However the par- 
ticular attention of the present work lies on the solution of the 
Bethe-Salper equation with the full off-shell dependence. Thus 
we will restrict the interaction kernel to a sum of contact terms, 
but refrain from the approximations mentioned above. 

Aside from the Weinberg-Tomozawa term, we will take into 
account the full set of meson-baryon vertices from the second- 
order chiral Lagrangian. The pertinent Lagrangian density was 
first constructed in 1 13] and reads in its minimal form Iil4i1 

+ bi/2(B[u^, W, B]^]) + bi{B[u^, {u", B}}} + b^iBBXu^u") 

+ !^5/6<ficr^''[[M/„ My], b]j + ibi{Bcj*"'u^){uyB) 

+ [«v, [D\B]]^]) + <B/[d„ [u\ [u„B]]^])) 

pl[(Br{u„ {«v, [D\ B]]}) + (Br[Dy, {u\ [u„ B]]])) 



2mo 

p^(2(By>'[D,,B]){u,u^} 



■ {ByfBXlDy, u^y + u^[Dy, u"])) 



(2) 



with the bi the pertinent dimension-two low energy constants 
(LECs). The LECs bo,D,F are the so-called symmetry breakers, 
while the bi (i - 1, . . . , 11) are referred to as dynamical LECs. 
On the one hand such terms may lead to sizable corrections 
to the leading-order result, see e.g. for the calculation of 
meson-baryon scattering lengths up to the third chiral order. On 
the other hand, including such terms with full off-shell depen- 
dence we hope to account for some of the structures created by 
the missing Born graphs. Let us denote the in- and out-going 
meson momenta by qi and q2, respectively. The overall four- 
momentum is given hy p = qi + p\ - q2 + Pi, where pi and p2 
are the momenta of in- and out-going baryon, respectively. Sep- 
arating the momentum space from the channel space structures, 
the chiral potential considered here takes the form: 

yii2-ivP) ^AwAqji +qb.) 

-H Ai4(^i ■ q2) + Asj[qi,q2] + Am 

+ Aui(^cf,{qi ■ p) + ^i(^2 ■ P)), 

where the first matrix Awt only depends on the meson decay 



constants F„_k, whereas A14, A57, A^n and Am also contain the 
NLO LECs as specified in appendix A. In going from the La- 
grangian Eq. (|2]i to the above vertex rule, we have left out some 
terms, which are formally of third chiral order The channel 
space is defined in accordance with the quantum numbers as 
well as the energy range of interest. For the purpose of gain- 
ing some insight on the nature of A(1405) it is spanned by six 
vectors corresponding to the following meson-baryon states: 
[K-p; K^n; 7r"A; ;r*'l*'; n+lT; 7r-E+). The influence of other 
much heavier channels will be absorbed in the LECs to be fitted. 

The strict perturbative chiral expansion is only applicable at 
low energies and certainly fails in the vicinity of (subthreshold) 
resonances. We extend the range of applicability by means of 
a coupled channel Bethe-Salpeter equation (BSE). Introduced 
in Ref. III6I1 it has been proven to be very useful both in the 
purely mesonic and also in the meson-baryon sector In con- 
trast to perturbative calculations this approach implements two- 
body unitarity exactly and in principle allows to generate res- 
onances dynamically. For the meson-baryon scattering ampli- 
tude r(^, ,fi',p) and the chiral potential V(^2' li ' P^ integral 
equation to be solved reads 

TitvlvP) ^V{q^,q^-p)+ 

\-^V{q^J-p)S{il>-mi)T{l,ci^-p\ 



■1 



(2nY 



(3) 



where S and A represent the baryon (of mass m) and the meson 
(of mass M) propagator, respectively, and are given by iS if) - 
iKp - m + ie) and iA(k) = ijik^ - + ie). Moreover, T, V, S 
and A in the last expression are matrices in the channel space. 

To treat the loop diagrams appearing in the BSE Eq. O), we 
utilize dimensional regularization. The purely baryonic inte- 
grals are set to zero from the beginning. In the spirit of our 
previous work IH, we apply the usual MS subtraction scheme, 
keeping in mind that the modified loop integrals are still scale- 
dependent. The scale fi reflects the influence of the higher-order 
terms not included in our potential. It is used as a fitting param- 
eter of our approach. 

The solution of the BSE Eq. (O with full off-shell depen- 
dence is obtained following the construction principles de- 
scribed in Ref. JSQ. As an extension of this approach we 
wish also to address an other issue here, namely analyticity. 
Let us first start with the one-meson-one-baryon loop-function 
^MBis = p^) in four dimensions. 



Imb(s) 



1 



1 



)4 f _M^(p- 1)2 



Applying the Cutkosky rules, one immediately obtains the 
imaginary part of this integral, given by -(^cms(/'^))/(87r V*) 
for q^^^ = {(s -(m + Mf)(s - (m - M)^)) /(4s). Keeping in 
mind the high energy behavior of the function Imb(s), we can 
obtain the real part of it via a subtracted dispersion relation 



Re(/MB(*)) = Re(/MB(*o)) ■ 



(s - So) p 



ds' 



, ImjlMBjs')) 
(s - s')(s- Sq)' 



2 



where Sthr - {M+nif and sq is a subtraction point chosen to not 
lie on the integration contour. The solution of the BSE Eq. (O 
corresponds to a bubble sum, containing exactly the same one 
loop-functions Imb- Thus one should in principle be able to 
write an equation similar to the last one for the scattering am- 
plitude 



Re(r(i)) = Re(r(io)) + 



Cds' 



Im(r(50) 



(4) 



where for the moment we have suppressed the ^ dependence. 
To the best of our knowledge, it is not possible to implement the 
last relation Eq. ^ into the BSE Eq. (O directly. To put it in 
other words, the BSE ansatz is known to produce poles on the 
physical Riemann sheet, which are forbidden by the postulate of 
maximal analyticity. Thus a scattering amplitude, which solves 
the Eq. does not satisfy Eq. (|4|. 

Nevertheless it is possible to find a solution of the BSE 
Eq. (O, which fulfills Eq. ^ at least approximately, as we wish 
to describe now. One way to do so is to keep only those solu- 
tions of the BSE, which do not produce poles on the first Rie- 
mann sheet 'near' the real (physical) axis. E.g. in Ref. jstl 
solutions producing poles for Im(Wcms - V^) < 250 MeV were 
excluded by hand. To overcome such unsatisfactory interven- 
tion into the fitting procedure we proceed differently. First, for 
a fixed s and sq we define the following quantity 



Re(T(s)-T(so))- 



Js, 



^^,Jm(7V))_^ 

S.hr (S-S')(S-S0) 



Rs(T(s)-T(so)) 



(5) 



Then the fitting parameters of our model are adjusted to min- 
imize the quantity A'full ~ ^msp /^'data' where the latter is 
based on the experimental data. It should be clear that such a 
procedure is not suited to overcome the unphysical poles. It en- 
sures, however, that they are moved far away from the real axis 
in a systematic manner, without manual intervention. This we 
consider an improvement of the model. 



3. Fit strategy 

We are now able to confront our approach with the experi- 
mental results. Throughout the present work we use the follow- 
ing numerical values (in GeV) for the masses and the meson 
decay constants: F„ = 0.0924, Fk = 0.113, = 0.135, 
M„H- = 0.1396, Mk- = 0.4937, M^o = 0.4977, lUp = 0.9383, 
m„ = 0.9396, ma = 1.1157, m^o = 1.1926, = 1.1894 
and niz- - 1.1975. The baryon mass in the chiral limit, mo in 
Eq. (|2]i, can be fixed to 1 GeV without loss of generality, as any 
other value only amounts to a rescaling of the unknown LECs. 

Secondly, for the experimental data we con- 
sider total cross sections for the processes K~p — > 
{K-p, K°n, n°lP, jt+S", ;r-2+) taken from 
Refs. iHElIillll. Moreover, we consider the following 



decay ratios 

r(K-p TT^I,-) 



7 = 
R„ = 
R, = 



= 2.38 + 0.04, 



TiK p 7r-S+) 

riK p 7r°A) 
T{K p — > neutral states) 

T{K p — > inelastic channels) 



0.189 ±0.015, 



= 0.664 ±0.011, 



where the first one is taken from Ref. 112 111 and the last two from 
Ref. ll22ll . Additionally to these quite old date we use a recent 
determination of the energy shift and width of the kaonic hydro- 
gen in the Is state, i.e. AE-iF/l = (283 ±42)-/(271 ±55)eV 
from the SIDDHARTA experiment at DAONE These are 
related to the K^p scattering length via the modified Deser-type 
relation 123 1 



AE — iT 12 — —la^fi^aK-p [l — latc paficQna — l)j , 

where a ^ 1/137 is the fine-structure constant, yUe is the reduced 
mass and a^-p the scattering length of the K p system. 

There are 17 free parameters in the present approach. First 
of all, the low-energy constants represent the heavy degrees of 
freedom of QCD, which are integrated out. Thus they have to 
be fixed in a fit to the experimental data. As a matter of fact, 
the fitting parameters of our approach correspond to the SU(3) 
low-energy constants, renormalized by the effects of the not 
included channels {?/A; 772"; A'+S ; /T^S"). Additionally, three 
subtraction constants have to be determined from a fit, which 
correspond to the logarithms of the undetermined regulariza- 
tion scales {/u^a, I^kn, M;rs)- 

To reproduce the experimental data as well as to preserve the 
property of analyticity as described in above, we minimize the 
following quantity ^^yj^j^ = Xmsp + Xqata^ where the first part 
is given in the Eq. (|5]l and the second part by the quantity 



2 ._ X 

^DATA 



d.o.f 



mi^ini- p^ ^ nC 



Here p is the numbers of the free parameters, n, is the num- 
ber of data points available for the observable / and is the 
number of observables. The present choice of ^^^^^ is crucial 
to ensure the equal weight of diff'erent observables, indepen- 
dently of the corresponding number of data points. The mini- 
mization itself is performed using MINUIT [24], especially the 
migrad strategy in two steps, which is due to the quite compli- 
cated structure of the BSE solution with full off'-shell depen- 
dence. First, parameters are found to minimize the ;^'pjjj^j^ in the 
on-shell parametrization. In the second step we turn on the "off"- 
shellness" slowly, minimizing in each step the;^'^^^^ and taking 
the parameters of the best fit from the previous step as starting 
values. Such a procedure guarantees preservation of the right 
analytic properties of the solution, found in the first step. 
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Figure 1: Total cross sections for the scattering of K^p to various channels versus the laboratory momentum. The black points with error bars denote the 
experimental data from 1 17, lR [T^[20h considered for the fits. Shaded (green) bands denote the lo" error bands calculated as described in the text. The reaction 
p — ► Aff" is not a part of our fit and presented here only for completeness. 



4. Results 

For the best fit x\,ata ~ 0-524 we obtain the following pa- 
rameter set (all bi in GeV"' and //, in GeV) 



log(yU/rA,/(lGeV)) = +1.155±0.181 
logOu;,5;/(lGeV)) = -0.008 ± 0.002 





log(//,A/(lGeV)) 




-0.010 + 0.003 




= +0.582 + 0.052 




= -0.332 ± 0.045 


bi 


= -0.310 + 0.092 


b9 


= +0.298 ± 0.087 


b3 


= +0.227 + 0.038 


bio 


= +0.198 + 0.058 


b4 


= -0.939 + 0.069 


bu 


= +0.516 + 0.058 


bs 


= +0.023 + 0.007 


bo 


= +0.710 + 0.211 


bt 


= +0.001 + 0.001 


bo 


= -0.291 +0.068 


bi 


= -2.518 + 0.110 


bp 


= -0.057 + 0.014 



We note that the LECs are all of natural size, indicating that 
all relevant physical mechanisms are included in the calcula- 



tions. The experimental data on total cross sections is repro- 
duced quite nicely, see FiglU We first wish to remark that due 
to insufficient number of data points the channel K p An^ 
is not considered as experimental input in the fit procedure. 
For completeness, we present the outcome of our approach for 
this channel in Fig[T] Secondly, all of the cross sections pre- 
sented here are due to the strong interaction only. Additionally, 
Coulomb interaction was taken into account in JSl |7[] via a non- 
relativistic quantum mechanical formula. Since this alone can- 
not count for an interference between the strong and the elec- 
tromagnetic interactions, we relegate the proper inclusion of the 
electromagnetic contributions to a future publication. 

The confidence bands presented in the above figure and the 
errors on further observables are calculated as follows: first we 
generate a large number (~ 10,000) of randomly distributed 
parameter sets in the error region given above. Then for each of 
these parameter sets we calculate the;t'Q^.j.^ and keep only those 
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Figure 2: Energy shift and widtli of kaonic iiydrogen as determined from the 
DEAR (51, the KEK [2Q and the SIDDHARTA [IJ experiments. The shaded 
area denotes the Icr region of our approach around the best fit value. 



sets, for which ;t'p^.p^ ~/?data - ^-^^^ Quantities calculated for 
these parameter sets are assumed to lie in the Icr region around 
the central value. 

The results for the threshold quantities are in excellent agree- 
ment with experimental data and read 

^E - ir/2 = +296!^^ - i 300!^^ eV, 



■V- +2 
R„ = +0.268!o;o8g . 



Rc = +0.643 



+0.015 
-0.019 



As a matter of fact, the shape of the 1 cr region for the energy 
shift and width of kaonic hydrogen cannot be assumed to be 
rectangular, see Fig. |2] The resulting scattering lengths for 
isospin 7 = and / = 1, i.e. aq and oi, are displayed in Fig. [3] 
in comparison to some older determinations and the determi- 
nation based on scattering data alone |5]. The inclusion of the 
SIDDHARTA data leads to much smaller errors, especially for 
fli. Our values for the scattering lengths are 



flo = -l-8l!|Si'^ + /0.92!°:^^fm, 
fli = -h0.48!J;:|? + / 0.87!°:^f; fm . 



(6) 



The inclusion of the Att" data in the fitting procedure could 
yield an additional constraint on the isospin 7=1 amplitudes 
and fix the value of ai as done in [Sd. We have not con- 
sidered this channel as an experimental input for the reasons 
given above. The scattering length for the elastic K p chan- 
nel reads ok-,, = -0.68![}[^ + / 0.90![J |^ fm. For compari- 
son, taking the SIDDHARTA data only, one obtains a^-p = 
-0.65!° j^ + i 0.81!° j|| fm, while Ikeda et al. fincfl OK-p = 
-0.70!°- + i 0.89!° fm. Therefore, these fundamental chiral 
SU(3) parameters can now be considered to be determined with 
about an accuracy of ~ 15%, 



Figure 3: Real and imaginary part of isospin and 1 KN — > KN scattering 
lengths. The light shaded (green) areas correspond to the Icr region of our 
approach around the central value (full circles). The darker (blue) areas corre- 
spond to the Icr region around central value (empty circle) from Ref. d. the 
cross and empty triangles denote older experimental values from |^ and iSHl . 
respectively. 



Having fixed the parameters of our model, we can extrapo- 
late the amplitudes of elastic K p scattering to the subthresh- 
old region, i.e. center-of-mass energies 1330 MeV < Wcms ^ 
1450 Me V. The result is presented in Fig.|4] For both real and 
imaginary parts of the amplitude the maximum lies close to the 
KN threshold and is quite narrow, which indicates the presence 
of a close-by pole. It is also worth mentioning that the error 
band gets smaller to low energies, different to the recent analy- 
sis by Ikeda et al. 

To obtain a more complete picture about the structure of 
A(1405), the amplitudes are analytically continued to the com- 
plex Wcms plane. Microcausality forbids poles on the first Rie- 
mann sheet, that is for Im(Wcm.s) > 0. This is fulfilled in our 
model automatically due to the restoration of analyticity as de- 
scribed above. On the other hand some pole structure has to 
be responsible for the functional form of the scattering ampli- 
tudes, see Fig.m Two poles are found on the second Riemann 
sheet for isospin 7 = 0, which is achieved via analytic contin- 
uation to Im(Wcms) < 0. We denote the second Riemann sheet 
connected to the physical axis in the region between the Ett and 
KN threshold as "R^jr and the one connected to the physical axis 
for Wcms > {Mfc + m^) as Hkn- We find that two poles lie on 
different Riemann sheets, the pole position reads 



Wi = 1428!^ - / 8!^ MeV . 

- / 75!^ MeV . 



'Rkn ■■ Wz = 1467+ 



'Here, the error bars are extracted from Fig. 4 of 



The real part of the position of the first pole agrees quite well 
with determination from Refs. [511101161171]. Its imaginary part 
agrees roughly with the determination of Refs. |5, 10] and is 
significantly smaller than extracted by Ikeda et al. 161 17|]. For 
the second pole, the situation is different, its imaginary part is in 
agreement with Refs. fsll^Ql, but the real part is much larger 
We have investigated the origin of these observations qual- 
itatively. First, from the analysis of nN scattering in the same 
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1340 1360 1380 1400 1420 1440 * 1340 1360 1380 1400 1420 1440 



WcMs [MeV] WcMs [MeV] 

Figure 4: Real and imaginary part of tlie K^p — » K^p scattering amplitude. The shaded band indicates the uncertainty of the calculation. The data point at 
Wcms = l^K + inp is determined from the energy shift and width of kaonic hydrogen from the SIDDHARTA experiment. 



framework, see Ref. jsl], it is known that off-shell effects can ac- 
count for large modifications of the pole positions. Setting the 
tadpole integrals to zero, we obtain immediately the solution of 
the BSE in the on-shell factorization. Note that this solution is 
still different to the one by Ikeda et al. smce no s-wave 

projection is performed. We found that in the present case the 
off-shell effects do not alter the pole position drastically. More 
precisely, the imaginary part of the first pole decreases and the 
one of the second increases by about 10 MeV. The real parts of 
both poles do not change significantly. Secondly, we noticed 
much smaller values of the NLO LECs found by Ikeda et al. 
additionally to the fact that the LECs (/ - 5, .., 11) were ne- 
glected there due to the s-wave projection. To keep track of this 
we scale down our LECs continuously from the values found 
above to zero. Such a solution of the BSE is of course by no 
means physical since no further fitting to experimental data is 
done here. Qualitatively, however, we observe that both poles 
move (the second one by about 100 MeV) to lower values of 
Re(Wcms)- The conclusion to be drawn is that difference in pole 
positions extracted in our approach and the one by Ikeda et al. 
is due to the differences in the fit strategies. 

Acknowledgments 

We are grateful to Peter Bruns for his stimulating remarks 
and cooperation. We thank Michael Doring and Wolfram Weise 
for comments. This work is supported in part by the DFG 
(SFB/TR 16 "Subnuclear Sti-ucture of Matter") and by the EU 
HadronPhysics3 project "Study of Strongly Interacting Mat- 
ter". 



+ /7i (</"'' [^^■"'", [A', A"]]) + {A'^'^IA', [Aj'KA"]])) 

+ h{{A'HA^\ {A', rm + {a'Ua', {a^\ r ))))), 

^sf" = bsiA'^i-^^KnA"]} + bdA'HlA^KA'lA"}} 

+ bn{{A''^A^'^){A'A'') + {A!'^A'){rA''))), 

+ b, (</"'■ [^^■"'", [^U"]]) + {A'hA\ [Aj\A"]]}) 
+ bg (</"'" [^^'1", {A', A"]]) + </"f [^', {Aj\ A"]])) 
+ bw{{A'HAjK {A', A"}}} + {A'Ha', {Aj\ A"}}})), 

A'f'" = IboiiA'U^XlA^'A'm) 

+ bD{{A''H{Aj'K {M, A'}}, A"}} + </"''{{i', {M, ^■'"'")), A"}}) 
+ fef {M, A'}}, A"]) + </"'"[{^', {M, /I-'"' )), A"]))), 

where A denote the 3x3 channel matrices (e.g. (p - (f/A' for the 
physical meson fields) and the F, are the meson decay constants 
in the respective channel. Moreover, M is obtained from the 
quark mass matrix M via the Gell-Mann Oakes Renner rela- 
tions, and given in terms of the meson masses as follows, M = 
idiag(M2 , - Ml, + , Ml, - Ml, + Ml, , MJ, + MJ„ - M^) . 



Appendix A. Couplings 

For the channel indices {b, j; i, a} corresponding to the pro- 
cess (piBa — > (pjBt, the relevant coupling matrices read 

A'„f"^-^{A'V,A'],A"]}, 
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